Generalization of the iterative perturbation theory and metal-insulator transition in 

multi-orbital Hubbard bands 
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The iterative perturbation theory of the dynamical mean field theory is generalized to arbitrary 
electron occupation in case of multi-orbital Hubbard bands. We present numerical results of doubly 
degenerate E 9 bands in a simple cubic lattice. The spectrum shows the electron ionization and 
affinity levels of different electron occupations. For sufficiently large Coulomb integral, a gap opens 
in the spectrum at integer filling of electrons and the system becomes insulator. The present scheme 
is easy to combine with the LSDA electronic structure theory. 

PACS numbers: 71.30.+h,75.40.Mg 
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The first principle electronic structure theory has been 
developed much recently in the field of strongly corre- 
lated electron systems, e.g. the LSDA+U method, Ji] 
the GW-approximation g etc. . The LSDA+U method 
can include the spin, charge and orbital fluctuation in 
real space, and the GW approximation is formulated to 
treat the dynamical correlation of electrons. Even so, 
these are still far from the goal since low energy excita- 
tion of local charge fluctuation, which plays an important 
role near the metal-insulator transition, is not properly 
formulated in the framework of the LSDA+U method or 
the GW-approximation. The metal-insulator transition 
and the anomalous metallic phase is now at the stage of 
our understanding in a unified picture of both low and 
high energy excitations due to application of the dynami- 
cal mean field theory (DMFT) developed by Georges and 
Kotliar. §!§ 

On the other hand, the DMFT might not be com- 
pletely established. Several methods have been devel- 
oped in order to calculate the Green's function in the 
DMFT. The DMFT projects a system onto the single im- 
purity Anderson model and this treatment is exact in the 
limit of the infinite dimension where the off-site Coulomb 
interaction can be neglected. The combination with the 
GW approximation is being developed to include the off- 
site Coulomb interaction. Q It is also tried to extend the 
single impurity approximation to the cluster approxima- 
tion. The effects of multi-orbitals were discussed [|[ ||] 
and the case of electron filling up to 1 was shown to be 
essentially the same as the case of the non-degenerate or- 
bital. [|| The combination with the LDA or the LSDA is 
also desirable. |h| [0], [H| 

To calculate the Green's functions of the DMFT, one 
can use several computational schemes such as the quan- 
tum Monte Carlo simulation (QMC), the iterative per- 
turbation theory (IPT), the non-crossing approximation 
(NCA), the exact diagonalization (ED), the numerical 
renormalization group (NRG) etc.. || QMC, ED and 
NRG are exact and, because of the exact calculation, 
would be limited for wide variety of applications. Other 
methods, IPT and NCA, are approximate methods and 
might be easier to combine with other calculation in re- 



alistic complex materials. The aim of this paper is to 
generalize the IPT to multi-orbital bands and to apply it 
to the metal-insulator transition in the doubly degener- 
ate E g orbitals. 

We start with the Hubbard-type Hamiltonian: 
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In the following calculations, we would simplify as 
U 3 mm , — U and J 3 mm , = J = 0, but its generalization 
is just straightforward. The local Green's function is 
then defined, since we neglect the k-dependence of the 
self-energy E as usual in the DMFT, as 
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where h(k) is the transfer matrix in k-space, /i is the 
chemical potential, £(iw) is the self-energy, and 1 is the 
unit matrix. In case of a paramagnetic and orbitally de- 
generate system, which is the case in the present paper, 
the self-energy E mo . :m ' ./(iw) and the local Green's func- 
tion G mCT:TO ' CT ' (iu>) are diagonal with respect to orbitals 
and spins and independent on them. 

The local Green's function G(iu>) can be then in a form 

of 



G(iui) = \ioj + n — A(iu>) — E(iu;)] 



(3) 



where A(iu>) is the effective hybridization function. From 
the definition of the "effective" medium, one defines an 
Green's function of the effective medium as 



A(iw)]' 



The second order self-energy can be defined as 
E( 2 >(r) - -U 2 (N dcg - l)G°(r) 2 G°(-r) , 



(4) 
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where iVd cg is the degeneracy with respect to spins and 
orbitals (here, in case of E g orbitals, JVd og = 2x2 = 4). 

The IPT was developed by Kajueter and Kotliar 
for non-degenerate orbital. Here, we assume a similar 
form of the self-energy as 



£M = Un(N dcg 1) + ^Mr. , 

1 — B(i(jj)l,^> (iu>) 



(6) 



where rid = ^2 ma nj m(T = nN deg and n is actually the 
occupation par each orbital. The coefficients A and B 
should be determined under the following guidelines: (1) 
A is determined so that the self-energy E(zw) is correct 
in the limit of iuj — > oo. (2) B(iu>) is determined so that 
the self-energy E(zu;) is correct in the atomic limit, i.e. 
fj-s-oo. 

The exact form of the self-energy in the limit of icu — > 
oo is generally given as 
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where ii/ is the electron-electron interactions. The 
brackets [• • •] and {• • • } are the commutator and anti- 
commutator, respectively. In the present case, the limit- 
ing form of the self-energy is 



lim icj{E 

iuj — >oo 
T2 
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Here we adopted a decoupling approximation 
{n ma n m i a i) ~ nn for (ma) ^ {m'a 1 ). Using the 
second-order perturbation, one gets a similar result as 

. lim iu E<£ !ro , a M - C^dcg - l)n°(l - n°) (9) 

where n° is the occupation number for the effective 
medium Green's function G° and, together with p,, is 
determined so as to satisfy the Luttinger theorem. Then 
the final form of A should be 



A = lim 

iuj — >oo 
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The required form of B{iu>) may be 
B(iu;) = S^ 8 '*^)- 1 - [E at (iw) - t/n(7V dcg 
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here E^^,^ is the second order self-energy in the atomic 
limit. 

The self-energy in the atomic limit can be calculated 
exactly, or even can be written down in an analytic form, 
through the atomic Green's function for the 7Vd eg -fold 
atomic level: p"5[ 
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, where p n ' d is the probability of finding an atom in a state 
of n' d electrons present. Here the atomic levels may be 
given as E = e d + Un' d according to the electron occupa- 
tion n' d = 0, 1, 2, • • ■ and all configurations are degen- 
erated when the electron occupation is identical. If we 
would adopt more realistic Coulomb and exchange inte- 
gral parameters U 3 mm , and J 3 mm i rather than the present 
simplified ones, the atomic levels should be split and they 
depend on the total angular momentum. 

In order to get a rigorous expression of the self-energy, 
one can assume the limit of the inverse temperature f3 — > 
oo with respect to U . Then the occupation probabilities 
of each configuration, for r < n c i < r + 1 (r : integer 
number), are 



Pr = l-(n d -r) , 
Pr+i = (rid - r ) , 
p r i =0 : r 7^ ?', r + 1 



(13) 



The resultant expression of B(iu>) may be as follows; 
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where, for r < rid < T + 1 (r : integer number), 

X = iuj — e' d — U r + fi , q = rid — r , 
TT 

a + b 



iuj — e d — Ur + fi 
° -{-q(N dcg -l)+r} , 
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Then finally we get an analytic expression of the self- 
energy of the present generalized IPT. The resultant ex- 
pression of the self-energy Eq. (@) with Eqs. ( 1C| ) and ( fl4| ) 
are actually simple extension of the IPT (lj and it is re- 
duced to the known expression when iVdcg — > 2. It is also 
similar to the results given in 10, ll) but not identical. 

We now present our results of the doubly degenerate 
Eg orbitals (A^og = 4). We consider the Slater-Koster 
type tight-binding Hamiltonian of the doubly degenerate 
E s orbitals in a simple cubic lattice with e d = 0. The ef- 
fective transfer integrals are assumed only between near- 
est neighbor pairs and V d d<x = 1/3, V d d^ = -2/3V d dcr, 
V d ds = l/GVdda- The relation among V dda , Vdd-K and 
VddS is the scaling properties of bare two-center inte- 
grals in the LMTO method. Here the band width is 
Q(V d da + Vdds) = 7/3 (the half band width D = 7/6). In 
a simple cubic lattice, there is no ddir interaction within 
the E ff orbitals and Vdd-K does not appear. Figure [l] shows 
the w(k)-curve and the real and imaginary parts of the 
local Green's function with E = and (3 — » oo. 

The transfer matrix {h(k)} is diagonalized at ev- 
ery k-point. The local Green's functions satisfy the 
relation G x 2 _ y 2 x 2_ y 2 {ilo) = G 3z 2_ r 2 t3z 2_ r 2(iuj) and 
G x 2_y2 3z 2_ r 2(iuj) = after the k-integration within the 
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FIG. 1: The energy bands w(k) along high symmetry lines 
(left) and the local Green's function (right) of Eg-bands in 
a simple cubic lattice. The high symmetry k-points are M 
(1,1,0), X (§,0,0), r (0,0,0), R (f , f , f ). The center of 
the band e d = and the transfer parameters are Vddu = 1/3, 
VddTT = — 2/Wdda, and VddS = Vdd^/d- The top E T and the 
bottom Eb of the bands are E T / B = e d ± 3( Vdda + Vdds)- The 
half band width is D = 3(Vdd<? + V dd s) = 7/6. 



whole Brillouin zone. To calculate Green's functions, 
8 x 343 k-points are used in the whole Brillouin zone 
and the generalized tetrahedron method for an arbitrary 
complex function is applied, which is the linear extrapo- 
lation scheme within a unit tetrahedron consisted of four 
k-points. The Pade approximation Jl7| is performed for 
analytic continuation of the Green's function from on the 
Matsubara frequency ioo n to the real w-axis. 

The imaginary part of the local Green's function is 
shown in Fig. |, for U = 1.0, 2.2, 3.0 and /3 = 30 with 
electron occupation < n < 1/2 or < rid < 2, where 
the energy region of lo < is occupied and that of u> > 
is unoccupied. The chemical potential decreases from 
l/2(7(iVdcg — 1) with decreasing n from 1/2. 



In case of n = 1/2 or rid = 2 (half-filled), the spectrum 
is very similar to that of the case of the non-degenerate 
band. The spectrum consists of the upper and lower Hub- 
bard bands with the electron-hole symmetry. At the crit- 
ical region (here, U = 2.2) one observes a sharp coherent 
peak at lo ~ 0. The change of the chemical potential with 
occupation is very rapid near n = 1/2 and n = 1/4 when 
U = 1.0 and U = 2.2 and it even jumps at these two occu- 
pations when U = 3.0. Therefore, the system becomes in- 
sulator at n = 1/4 and 1/2 when U = 3.0 (U/D ~ 2.57). 
These occupations correspond to the cases of the total 
electron numbers rid = 1 and 2, respectively, and this 
is actually the case where some Hubbard sub-bands are 
completely filled and others are empty. The band gap of 
U = 3.0 at n = 1/4 is narrower than that at n = 1/2 and 
then the critical U / D-value for metal-insulator transition 
may be smaller at n = 1/4 than that at n = 1/2. This 
prediction is consistent with the results by QMC and, 
furthermore, it is presumed that the critical U / D-vahie 
is smaller than the model of semicircular density of states 



in a multi-orbital case ( (U/D) cr = 3.1 at n — 1/2). || 
This fact suggests that the critical U/D- value would de- 
pend upon the Lo(k) structure or the shape of the density 
of states. In insulating case, U — 3.0 with n = 1/4 
and n = 1/2, the self-energy shows the singular behavior 
£ ~ J— near iuj n ~ 0. On the contrary, in almost all 
cases, tfie system is metallic and the self-energy behaves 
as S — Un(Ndcg — 1) ~ iio n which is the property of the 
Fermi liquid. The case U = 2.2 is at the critical region of 
the metal-insulator transition, and one can see very sharp 
coherent peak at lo = when n = 1/4 and n = 1/2. In 
this region, the self-energy shows the common feature of 
the Fermi liquid. 

With small hole doping, the spectrum changes very 
drastically especially when U — 2.2 and U = 3.0. A 
sharp peak appears at U below the lower Hubbard band, 
whose width depends upon the value U or the separation 
from the main lower Hubbard band. In this region of 
1/4 < n < 1/2 or 1 < rid < 2, the lower energy state is 
the mixture of the ground states of rid = 1 and rid = 2. 
Then, the ionization and affinity levels for both rid = 1 
and 2 configurations should appear in the spectrum. The 
satellite structure of all configurations iid = 0, 1, 2, ••• 
could appear, but the intensity of satellites would depend 
on the energy separation. These satellite structure orig- 
inates from the self-energy correctly given in the atomic 
limit. 

Near the occupation n ~ 1/4, the satellite at U above 
the upper Hubbard band shrinks. When the concentra- 
tion n crossed the value 1/4, i.e. < rid < 1, the upper 
satellite disappears and the lower satellite merges into the 
main (lower) Hubbard band, since the main contribution 
to the spectrum comes from the mixture of rid — and 
1, rather than rid = 1 and 2 of the case 1 < rid < 2. In 
fact at this concentration, the chemical potential shifts 
rapidly or jumps and the satellite grows to a main lower 
Hubbard band. Then the behavior for < n < 1/4 
is rather similar to the case of < n < 1 of the non- 
degenerate orbital and the difference is only the fact that 
the intensity of the lower and upper bands starts from the 
ratio of 1 : 3 near n ~ 1/4 rather than l:lofn.~l/2 
in case of the non-degenerate orbit. 

The k-dependent spectrum can be analyzed through 
the k-dependent Green's function G mm ' a {Lo, k) = [(lo + 
fi) 1 - e° d 1 - h (k) - S ( lo ) } o . The band dispersion width 
maybe sensitively depend upon the electron occupation. 
In the region of w ~ 0, the Green's function behaves 
as G(w,k) ~ Z/{lo1 — Zh(k)}, where Z is the renor- 
malization factor. This situation may be seen, in case 
of the small doping from half filling, the top of the co- 
herent peak shows a small dip which is a characteristic 
feature of the lattice Green's function of the E 5 state. 
The height of the coherent peak is almost equal to each 
other at n ~ 1/2 and n ~ 1/4. Then the renormalization 
factor or the effective mass may be also equal. 

In general, it may be very important to have a sim- 
ple scheme for solving the self-consistent Green's func- 
tion in multi-orbital cases, since further generalization is 
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FIG. 2: The electron occupation dependence of the imaginary part of the local Green's functions for the doubly degenerate 
E 3 band of the half band width D = 7/6. (a) U = 1.0, (b) U = 2.2 and (c) U = 3.0. The numbers in the figure denotes the 
electron occupation n for each spectrum. The inverse temperature /3 = 30. The region co < is occupied and that to > is 
unoccupied. In case of the electron occupation n = 1/4 and 1/2 at U = 3.0, a gap opens and systems are in insulating phase. 



needed for application to realistic materials showing very 
interesting interplay of spin, charge and orbital ordering. 
For that purpose, the present scheme is very practical 
and straightforward and, furthermore, the generalization 
for magnetically ordered systems and clusters containing 
several atoms may be quite simple. 

In conclusion, we developed a new generalization of 
the self-energy of the DMFT-IPT applicable to arbitrary 
electron occupation 0<n<lor0<n d < N dcg . 
The spectrum shows the electron ionization and affinity 



levels of different electron occupations. For sufficiently 
large Coulomb integral, the system becomes the insulat- 
ing state at integer filling of n^. This generalization of 
the IPT may be easily applicable to more general cases 
and the combination with the LSDA in more complex 
systems. 
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